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Abstract 

This paper gives an integrable hierarchy of nonUnear evolution equations. 
In this hierarchy there are the following representative equations: 

Ut = 

"'t fJx ! 

u 

'^xxt ~t~ '^^xx^x ~t~ ^xxx'^ 0* 

The first two are in the positive order hierarchy while the 3rd one is in the neg- 
ative order hierarchy. The whole hierarchy is shown integrable through solving 
a key 3x3 matrix equation. The 3x3 Lax pairs and their adjoint representa- 
tions are nonlinearized to be two Liouville-integrable canonical Hamiltonian 
systems. Based on the integrability of 6A^-dimensional systems we give the 
parametric solution of the positive hierarchy. In particular, we obtain the 
parametric solution of the equation ut = d^u~3 . Moreover, we give the trav- 
eling wave solution (TWS) of the above three equations. The TWSs of the 
first two equations have singularity and look like cusp (cusp-like), but the 
TWS of the 3rd one is continuous. For the 5th-order equation, its parametric 
solution can not include its singular TWS. We also analyse the Gaussian ini- 

_ _ 2 

tial solutions for the equations ut = d^u s , and Uxxt + -iUxxUx + UxxxU = 0. 
One is stable, the other not. Finally, we extend the equation ut = d^u"^ to 
a large class of equations ut = c)^u~"^/", Z > 1, n 7^ 0, 111,71 £ 7^, which still 
have the singular cusp- like traveling wave solutions. 

Keywords Hamiltonian system, Matrix equation, Zero curvature representation. 
Parametric solution. Traveling wave solution. 
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1 Introduction 

The inverse scattering transformation (1ST) method plays a very important role 
in the investigation of integrable nonlinear evolution equations (NLEEs) 0. This 
method has been successfully applied to solve the integrable NLEEs in the form of 
soliton solutions. These NLEEs include the well-known KdV equation |r5|, which 



is related to a 2nd order operator (i.e. Hill operator) spectral problem |T6|, |T8|, the 
remarkable AKNS equations 0, which is associated with the Zakharov-Shabat 
(ZS) spectral problem and other higher dimensional integrable equations. 



In the theory of integrable system, it is significant for us to search for as many 



new integrable evolution equations as possible. Kaup [|12| studied the inverse scat- 
tering problem for cubic eigenvalue equations of the form ip^xx + QQi^x + 6Rip = Xip, 
and showed a 5th order partial differential equation (PDE) Qt + Qxxxxx + 'iO{QxxxQ + 
\QxxQx) + ISOQsQ^ = (called the KK equation) integrable. Afterwards, Kuper- 
schmidt fl^ constructed a super-KdV equation and presented the integrability of 
the equation through giving bi-Hamiltonian property and Lax form. Very recently, 
Degasperis and Procesi proposed a new integrable equation (called DP equation) 
with the soliton solution of peakon type. The DP equation is an extension of the 
Camassa-Holm (CH) equation @], and is proven to be associated with a 3rd order 
spectral problem 0: ''\)xxx = '^x~ Am^^ and to have some relationship to a canonical 
Hamiltonian system under a new nonlinear Poisson bracket (called Peakon Bracket) 



IC]. In [11], we extended the DP equation to an integrable hierarchy and deal with 
its parametric solution and peaked stationary solutions. 

In Ref. the authors studied the DP equation, of a similar form to the 
Camassa-Holm shallow water wave equation, and proved the exact integrability 
of this equation by constructing its Lax pair. The DP equation is related to a 
negative flow in the Kaup-Kupershmidt hierarchy via a reciprocal transformation. 
The infinite sequence of conserved quantities is derived together with a proposed 
bi-Hamiltonian structure. The equation admits exact solutions in the form of a 
superposition of multi-peakons, and is compared with the analogous results for 
Camassa-Holm peakons. 

The present work is motivated on the basis of the eigenvalue problem -02. — 
o?"^xxx = a'^\rm/j, which was introduced by Desgaperis, Holm and Hone 0. Here, 
we consider the limit case of a going to infinity. That is, we get the 3rd order 
spectral problem = —Xrml). In this paper, starting from that problem, we give 
an integrable hierarchy, and through solving a key matrix equation we explicitly 
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provide the Lax pairs for the whole hierarchy. The following equations 



Ut = d^u 



lu i. 



Ut = ' (2) 



u 



'^xxt + Sltjjajltj; + Uxxx'U — 0, (3) 

are three representatives in the hierarchy. The first equation is a reduction of some 



2 + 1 dimensional equation [jl3l . The second one is new, and the third one is actually 
twice derivatives in x of the Rieman shock equation 0. 

Konopelchenko and Dubrovsky [|13[ already pointed out equation (|1|) is integrable 



from the reduction view point [|TB[ , but did not discuss the spectral problem and the 
solution of the equation. Here we deal with its spectral problem and representation 
of solution from the constraint view point. We give the parametric solutions for the 
hierarchy, particularly for equation Furthermore, we obtain the traveling wave 
solution (TWS) for equations (|I|), (0), and (0). The first two look like a class of 
cusp sohton solutions (called 'cusp-like', but not cusp soliton |2^). The TWSs of 
equations (0) and @ have singularity, but the TWS of equation (H) is continuous. 
Additionally, for the 5th-order equation (Q), its smooth parametric solution can 
not include its singular TWS. Equation has the compacton-like and parabolic 
cylinder solutions. We also analyse the Gaussian initial solutions for equations 
Ut = d^u~3 and u^^t + '^u^^u^ + u^xxU = 0. The former is stable, and the latter not 
(see Figures 7, 1 - 5). Finally, we extend the equation ut = d^u~3 to a large class 
of equations Ut = I > 1, n ^ 0, m,n E "L, which still have the singular 

cusp-like traveling wave solutions. 

The whole paper is organized as follows. Next section is saying how to connect 
the above three equations to a spectral problem and how to cast them into a new 
hierarchy of NLEEs. In section 3, we construct the zero curvature representations 
for this new hierarchy through solving a key 3x3 matrix equation. In particular, 
we obtain the Lax pair of equations (|l]), (0), (|^), and therefore they are integrable. 
In section 4, we show that the 3rd order spectral problem related to the above three 
equations is nonlinearized as a completely integrable Hamiltonian system under 
some constraint in M^^. In section 5 we give the parametric solution for the positive 
order hierarchy of NLEEs. We particularly get the parametric solution of equation 
(|1]). Moreover, in section 6 we obtain the traveling wave solutions for equations (|I|), 
(^, and (^, and we also analyse the Gaussian initial solutions for the equations 
Ut = d^u~s and u^xt + 'iuxx^x + UxxxU = 0. We propose a family of new equations 
Ut = / > 1, n 7^ 0, m,n E Tj, which still have the singular cusp-like 

traveling wave solutions. Finally, in section 7 we give some conclusions. 
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2 Spectral problems and a new hierarchy 

Let us consider the following 3rd order spectral problem 



i>xxx = -Xwi/j (4) 

and its adjoint problem 

i^*xxx = Xuij*. (5) 

Take u ^ u + e5u, A ^ A + e5X, and ^|e=o- Then, we have their functional 
gradient |^ with respect to the potential u 

5X XiJtp* _ VA 



6u E E 

where 



(6) 



VA = Xi:^*, (7) 
E = uipip*dx = constant, 



and Q = (— cxd,oo) or f2 = (0,T). In this procedure, we need the boundary condi- 
tions of u decaying at infinities or of u being periodic with period T. Usually, we 
compute the functional gradient |^ of the eigenvalue A with respect to the potential 
u by using the method in Refs. ^ 



Through doing five times derivatives of Eq. (|^), we find 

(VA)™ = -3X\2ud + du)i^rx-r^x), 
[iljipl - ij*^px)xxx = {ud + 2du)VX, 

which directly lead to 

KVX = AVVA, (8) 

where 

K = d\ (9) 
J = -3{2ud + du)d-%ud + 2du). (10) 

Hint: Here we do not care about the Hamiltonian properties of the operators 
K, J, but need 

K-^ = d-\ 

= -lt,-2/3a-l^-V353^-l/35-l^-2/3_ 

27 
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They yield 

£ = J-^K = -1^,-2/39-1^-1/353^-1/3^-1^-2/3^5 (11) 
27 ' V ; 

£-1 = R-^j = -3d~^2ud + du)d^^{ud + 2du). (12) 

By this pair of operators we define the hierarchy of nonhnear evolution equations 
associated with the spectral problems (^ and (^). Let Gq G Ker J = {G & 
C°°(M) \ JG = 0} and G-i E Ker K = {G e C°°(M) | fsTG = 0}. We define the 
Lenard sequence 

j<0, jez. ^^'^ 

where £ = J~^K is called the recursion operator. Therefore we produce a new 
hierarchy of nonlinear evolution equations (NLEEs): 

ut, = JGk, e Z. (14) 

Apparently, this hierarchy includes the positive order [k > 0) and the negative 
order {k < 0) cases. Let us now give several representative equations in the hierarchy 

• Choosing G^i = | G Ker K yields the first equation in the negative hierarchy: 

ut + vu^ + "iv^u = 0, M = v^^. (15) 

This equation is actually: v^xt + 'iv^x^x + v^xxV = which is equivalent to 
d^ivt + Wx) = 0. Obviously, v = cix + cq (ci, cq are two constants) is a 
special solution of this equation. In section 6, we will study its traveling wave 
solution. 

The second equation in the negative hierarchy is: 

Ut_2 + 'iiuxW + 3uWx) = 0, 

Wxxx + -{2U^V + 3uV:r) = 0, 

2 

• Choosing Gq = u~3 g Ker J leads to the second equation in the positive 
hierarchy: 

Ut = d^u'^. (16) 

Konopelchenko and Dubrovsky ever pointed out that this equation is inte- 
grable and is a reduction of some 2 + 1 dimensional equation [O. But they 



did not study the solution of this equation. In the following, we study the 
relation between this equation and finite-dimensional integrable system and 
will find that it has parametric solution as well as the traveling wave solution 
which looks like a cusp. 

The third equation in the positive hierarchy is: 

Wx = U-2/3(^-2/3^)^^_ 

• Choosing another element Go = ^" ^)xx-2{u s)^ ^ j^^^ j gj^^g following 
representative equation in the positive hierarchy: 

u 

This equation also has a singular cusp-like traveling wave solution. 

Of course, we may produce further nonlinear equations by selecting other elements 
from the kernels of J, K. In the following, we will see that all equations in the 
hierarchy (0) are integrable. Particularly, the above three equations (|15|), 
(P!B|), (pr7|) are integrable. 

3 Zero curvature representations 

Letting -0 = -01, we change Eq. to a 3 x 3 matrix spectral problem 





f/(n, A) = 1 , ^ = . (19) 



Apparently, the Gateaux derivative matrix U^:{^) of the spectral matrix U in the 
direction ^ G C°°(M) at point u is 




U{u + ei)= \ I (20) 

£=0 



which is obviously an injective homomorphism, i.e. ?7*(.^) = -v^ = 0. 

For any given C°°-function G, we construct the following 3x3 matrix equation 
with respect to ^ = V{G) 

- [[/, V] = U,{KG - X^JG). (21) 
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Vi2 


^13 


V21 


V22 


^23 


V^31 


V^2 


V^33 



Theorem 1 For the spectral problem ( |7^ and an arbitrary -function G, the 
matrix equation i^2l\ ) has the following solution 

-G" - 3Xd-^TG 3(G" + Xd^^G) -6G 
V = X\ -G'" -3Xd-^uG' 2G" 3{-G' + Xd-^TG) |, (22) 

-G"" - SX^ud-^G G'" - 3Xd~^uG' -G" + ^Xd^^TG 

where d = dx = -§^-, T = ud + 2du, and the supscript ' means the derivative in x. 
Therefore, J = —3T*d~^T (T* is the conjugate off). 

Proof: Set 



V 



and subsitute this into Eq. (|2TD. This is a overdetermined equation. Using some 
calculation techniques ||19[, we obtain the following results: 

V^i = -XG" -3X^d'^TG, 

V12 = 3{XG' + X^d-^TG), 

Vis = -QXG, 

V21 = -XG'" -3X^d-\G', 

V22 = 2XG", 

V23 = 3X{-G' + Xd-^TG), 

V31 = -XG"" ~3X\d-^TG, 

V32 = XG'" -3X^d-^uG', 

V33 = -XG" + 3X^d~^TG, 



which completes the proof. 

Theorem 2 Let Gq G Ker J, G-i G Ker K, and let each Gj be given through Eq. 
g^. Then, 

1. each new vector field Xk = JGk, A; G Z satisfies the following commutator 
representation 

Vk,x-[U,Vk] = U,{Xk), VA;gZ; (23) 

2. the new hierarchy (\m), i.e. 

ut^ = Xk = JGk. \Jk G Z, (24) 
possesses the zero curvature representation 

Ut, - Vu^x + [f/, Vu] = 0, Vfc G Z, (25) 
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where 



' A; > 0, 

V, = Y.V{G,)\^^''~^~^\ Yl = { 0''"°' ^ = 0' (26) 



and V{Gj) is given by Eq. (WW with G = Gj. 
Proof: 

1. For A; = 0, it is obvious. For k < 0, we have 

- [U, V,] = - ^ - [U, V{G,)])X'^'-^-'^ 

j=k 

= -J2U, {KG, - X^KG.^i) a2(^--^'-i) 

j=k 

\j=k 

= U^KGk^i-KG^iX^') 
= U^KGk-i) 
= U^Xk). 

For the case of A; > 0, it is similar to prove. 

2. Noticing Ut^ = U^{ut^), we obtain 

Ut^-Vk,, + [U,Vk] = U,{ut,-Xu). 
The injectiveness of U^^ imphes item 2 holds. 
From Theorem ^, we immediately obtain the following corollary. 
Corollary 1 The new hierarchy ([7^ has Lax pair: 



i^xxx = -Xw/j, (27) 
= 5ZA2('=^^)-i[-6G,^/.,, + 3(G;. + Aa-=^TG',)V',-(G;' + 3Aa^=^TG',)^^^ 

(28) 

where the related symbols are the same as Thereom ^ and Thereom 



So, all equations in the hierarchy (|14|) have the Lax pair and are therefore inte- 
grable. In particular, we have the following special cases. 



When we choose G-i = |, equation (|T^) has the following Lax pair: 

= U{u,X)^, (29) 
= Viu,X)^, (30) 



where u = Vxx, U{u, A) is defined by Eq. (|T9D, and V{u, A) is given by 

(Vx —V A~^ \ 
-V \. (31) 
Xvu —Vx I 

Apparently, Lax pair (|29| ) and (0) is equivalent to 

'^xxx = -Xuifj, (32) 
i^t = A'Vxx - vijjx + Vx^. (33) 

which is a limit form in Ref. H when a goes to oo. 



In a similar way, choosing Gq = u i gives the Lax pair of equation (|T^), i.e. 

(^ ^xxxxx 

i^xxx = -Xuxjj, (34) 

2 2 2 

i{jt = -QXu'ii/jxx + SX{u~3)^iJj^ - X{u~3)xxip- (35) 



This Lax pair is different from/inequivalent to the result in Ref. |[13|| . 



Furthermore, through choosing Go = ^)xx ^2{u k)^ ^ ^j^^^ ^j^^ 

equation (^ has the Lax pair: 



i^xxx = -Xuip, (36) 
ijt = -6XGoiJxx + ^XiG'o + 3Xu"'^)iJjx-X{G'^ + 9X{u--^)xx)^- 

(37) 

4 Nonlinearized 6A^-dimensional integrable system 
from spectral problems 



To discuss the solution of the hierarchy (p!4D, we use the constrained method which 
leads finite dimensional integrable system to the PDEs (|14D. Becasue Eq. (|^) / (p!8|) 
is a 3rd order eigenvalue problem, we have to investigate itself together with its 
adjoint problem when we adopt the nonlinearized procedure 0. Ma and Strampp 
^ ever studied the AKNS and its its adjoint problem, a 2 x 2 case, by using the 



so-called symmetry constraint method. Now, we are discussing a 3 x 3 problem 
related to the hierarchy ([T^). 
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Let us return to the spectral problem (|I^) and consider its adjoint problem 
and change it to the following matrix form 




u\ 

= I -1 I**, = I -0; I , (3J 
,0-10 

where i/j* = ip^. 

Let Xj {j = 1, N) he N distinct spectral values of ( ]T8| ) and (|38|), and qij, q2j, qsj 
and pij , p2j , P3j be the corresponding spectral functions, respectively. Then we have 

qix = q2, 

q2x = qs, (39) 

q-sx = -uAqi] 

and 

P2x = -Pi, (40) 

P3x = -P2, 

where A = diag{Xi, X^), qu = (?fci, 9^2, ^fcTv)'^, Pk = {Pki,Pk2, ■■■,PkNV , k = 
1,2,3. 

Let us consider the above systems in the whole symplectic space (M^^, dpAdq). 
We directly impose the following constraint: 

N 

«-i = 5^VA„ (41) 
i=i 

where VAj = Xjqijp^j is the functional gradient of Xj for spectral problems (|T8D and 
Then Eq. (^) is saying 

u= {Aqi,p3)~'^ (42) 



which composes a constraint in the whole space M . Under this constraint, Eq. 
( p9D and its adjoint (^) are cast in a Hamiltonial canonical form in M^^: 

qx = {q,H+}, 

Px = {p,H+}, ^ ^ 



with the Hamiltonian 



= {q2,Pi) + {q3,P2) + ^=^=, (44) 

V(Agi,P3) 
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where p = {pi,p2,pzY , q = {qi^qiAzY ^ 



■, ■) stands for the standard inner 



product in 
follows: 



and we modify the usual Poisson bracket of two functions Fi, F2 as 



i=i ^ 



dqi dpi 



dFi dF2 
dpi ' dqi 



(45) 



which is still antisymmetric, bilinear and satifies the Jacobi identity. 

To see the integrability of the system (^31), we take into account of the time part 

= 14.* and its adjoint = -Vf where Vk is defined by Vk = Ejl^ V{Gj)X'^^''-^-^\ 

and V{Gj) is given by Eq. (H) with G = Gj. 

Let us first look at Vi case. Then the corresponding time part is: 



X 



( 



\ 



-{u §; 

-(U~3), 



2 N 

3 

--^ 

3 



3{u- 
2{u- 



—6u 3 \ 

-3(w"i)a: 



_2 ' 
[U 3 



(46) 



XX j 



and its adjoint part is: 



X 



\ 



-3(^-3), 
6u~3 



[u 3 



6 Am 3 



\u 3 



— [U 3 



XXX 

U 3 



QXu3 



(47) 



3(m 3)^ 

Noticing the following relations 

u^i = {Aqi,p3)~^ , 

(Ag3,P3) + (Agi,Pi)-2(Ag2,P2), 
3((Ag2,Pi) - (Ag3,P2)), 

6 {Aq3,pi) + 3 (Agi,p3)"^ ( (A'gi,P2> + (A'g2,P3> 

1), and cast the non- 



(48) 



(w-i; 

fw~3)^_^_^ 



we obtain the nonlinearizations of the time parts (|46|) and 
learized systems into canonical Hamiltonian system in M^^: 

qt,={q,F+}, 
Pti = {p,Fi}^ 



with the Hamiltonian 

= ( {AquPi) + (A53,|>3) + 2 (Ag2,P2) ( (Agi,Pi) + {Aq3,P3) - {H2.P2) 
+3((Ag2,P3) - (Agi,p2))((Ag2,Pi) - {Aq^^p^)^ -6{Aqi,p3) {Aqs^p,) 

+ rj^ , ( (A'gi,P2> + {A\2.Pz) ) . (49) 

A direct computation leads to the following theorem. 
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Theorem 3 



{H+,F+} = 0, (50) 

that is, two Hamiltonian flows commute in M^^. 

Furthermore, for general case V^, A; > 0, A; G Z, we consider the fohowing 
Hamihonian functions 



j=0 
k—1 

j=0 
k—l 

+3j2{{A'^^%,Ps) - {A'^'-'qi,P2)){{A'^'-^^'%,Pi) - {A'^'-^^-%,P2) 

j=0 
j=0 



j=0 
k 



-3 E ( (A^^g2,P3> + {A''qi,P2) ) ( {A'^'~'k2,Pl) + {A'^'~^k3,P2) 
j=0 

+3H+^{A''q„p,) + {A'%,Ps)). (51) 
Then through a lengthy calculation, we find 

F+} = 0, F+} = 0, fc, / = 1, 2, ... . (52) 

That is. 

Theorem 4 ^4// canonical Hamiltonian flows (Fj^) commute with the Hamiltnonian 
system (]^^. /n particular, the Hamiltonian systems anc? are compatible 
and therefore integrable in the Liouville sense. 

Remark 1 In the proof of this Theorem, we use the following two facts: {qi,P2) + 
{q2,P3) = ci, and {qi,pi) — {qs^Ps) = c^- They always hold along x—flow in the 
whole . Here ci , C2 are two constants. 



Remark 2 In fact, the involutive functions are generated from the nonlineariza- 
tion of the time part = Vk"^ and its adjoint part = —Vlf'^* under the con- 
straint where Vk is defined by Vk = Yl'jZo^{C!j)X'^^''~^~^\ and V{Gj) is given 
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by Eq. (^Wj with G = Gj. In this calculation process, we use the following equalities: 





G, 


= -(A2^+V,P3>, J = 0,1, 2... 








/a 21+1^ ^ \ /A2i+1„ ^ \ 








— \^»- 93,^3/ + ~ 


- ^ q2,P2/ 






of /a 21+1 „ „ \ /a 21+1^ „ 

= (A -^^ g2,Pi) - (A -^^ g3,p2 








= 6(A2^+ig3,Pi>+3(Agi,p3)"^ 




d- 


'mG'^ 


= (A2^g3,p2> + (A'^g2,Pi>, 






~'TG, 


= (A2^gi,pi>-(A2^-g3,P3>, 




d- 


-'TG, 


= -((A2^gi,p2> + (A'^g2,P3>). 





5 Parametric solution 

Since the Hamiltonian flows (H^) and {Ff!') are completely integrable in and 
their Poisson brackets {H^,Fi^} = {k = 1,2,...), their phase flows g^+ 

commute 0. Thus, we can define their compatible solution as follows: 

q{x,tk) \ _ u ( q{x^A) \ , _ ^ 2 r53) 
p{x^t,))-'^^'K\p{x\tl))^ fc-1,2,..., (53) 

where t° are the initial values of phase flows g^+, (7^+. 

Theorem 5 Let q{x,tk) = {qi, q2, qs)^ , p{x,tk) = {pi,P2,P3)'^ be the common solu- 
tion of the two commutable Hamiltonian flows {H+) and (F^) in MP^ . Then 

u = ^ (54) 

y {Aqi{x,tk),P3{x,tk)f 

satisfies the positive order hierarchy 

ut, = Jd'-u-l, A; = 1,2,..., (55) 

where the operators C = J~^K, J, K are given by Eqs. ( |7^ j and respectively. 

Proof: Direct computation completes this proof. 

Theorem 6 Let p{x,t),q{x,t) (p{x,t) = {pi,P2,P3V ,q{x,t) = {qi,q2-,qzY) be the 
common solution of the two integrable commutable flows and ^^ ), then 

u = (56) 
y {Aqi{x,t),p3{x,t)f 

satisfies the equation: 

ut = dlu-i. (57) 
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Proof: Doing five times derivatives in x on both sides of Eq. (|5B|), we obtain 
dlu-l = 9w((A2g3,p3) - (AV,pi)) +3M,.((A2gi,p2> + (A2g2,P3>), (58) 



where 



3 ( (AV,P2) + {A^q2,P3) ) ( {Aq2,P3) - (A?i,P2) 



2 (Agi,P3) 

On the other hand, doing the derivative in t on the both sides of Eq. (^) yields 

3 (Ap3,gi) + (Agi,p3) 
Ut = u n ; 

2 (Agi,p3) 

3 (A.3,ff)-(A.,f^) 
= — li— ^ ^ — 

2 (Agi,P3) 

Substituting the expression of into the above equahty and calculating, we find 
that this final result is the same as the right hand side of Eq. (|58| ) , which completes 
the proof. 

6 Traveling wave solutions 

First, Let us compute traveling wave solution for equation (^. Set u = f{^), ^ = 
X — ct {c is some constant speed), then after substituting this setting into equation 
we obtain 

-cf" + 3/7' + /'7 = o, 



I.e. 



if - 2c/)'" = 0. 



Therefore, 



(/ - cf = + B^ + C, y A, B, C eR. (59) 
So, the equation Uxxt + ^UxxUx + UxxxU = has the following traveling wave solution 

u{x, t) = c± ^/A{x - cty + B{x - ct) + C. (60) 
Let us discuss several special cases: 
• When c = 0, we get stationary solution 



U[X] 



±VAx2 + Bx + C, VA, 5, C e M, (61) 



which can be a straight line, circle, ellipse, parabola, and hyperbola according 
to different choices of constants A, B, C. 
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When c 7^ and A 7^ 0, then we have 



A[x-ct^^ + ^ , VA, 5, C e M, (62) 

therefore if 4AC — = this solution becomes 

u{x, t) = c^\fA[x- ct^ — j , Vy4 > 0, 5 e M (63) 
For example, setting c=l, ^4 = 1, i? = yields 

M(a;, t) = 1 — |x — t|, (64) 

and 

u{x,t) = \^\x-t\. (65) 



The former looks like a compacton solution pO|, ^. The latter is a "V"-type 
solution. 

When c 7^ and A = 0, then we have 

w(x, t) = c± ^B{x - ct) + C, y B, C eR, (66) 

which is a parabolic traveling wave solution if S 7^ and becomes a constant 
solution if i? = 0. For example, the following 



u{x,t) = l + y/x-t,x-t>0, (67) 

and 

u{x,t) = 1- Vx-t,x-t>0. (68) 
are two special solutions. 

So, the 3rd-order equation u^xt + '^UxxUx + UxxxU = has the continuous traveling 
wave solution (pOl). 

The Gaussian initial solution of this 3rd-order PDE is stable from t = to t = 64 
(see Figure p.But after t = 64 this solution is not satble (see Figures 2-5). 
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Figure 1: Stable period from t = to t = 64 for the equation Uxxt + '^UxxUx + UxxxU = 
under the Gaussian initial condition. This figure is very like the Burgers case 
Uxxt + "iUxxUx + UxxxU — ^Uxxxx = which is formed through adding small viscosity 
term euxxxx to the equation. For instance, when e = 0.01, the equation has Figure 



equation ujxxt) + 3u_(xx)u_x + u_{xxx)u = 



20 30 40 50 60 70 




90 100 



Figure 2: Solution shape at t = 70 for the equation Uxxt + '^UxxUx + UxxxU = 0. 



equation ujxxt} + 3u_{xx)u_x + u_{xxx}u = 




90 100 



Figure 3: Solution shape at t = 80 for the equation Uxxt + "iUxxUx + UxxxU = 0. 
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equation ujxxt} + 3u_(xx)u_x + u_{xxx}u = 



10 20 30 40 50 



70 80 90 100 



Figure 4: Solution shape at t = 82 for the equation Uxxt + 'iUxxUx + UxxxU = 0. 



equation ujxxl) + 3u_{xx)u_x + u_{xxx)u = 



ID 20 30 40 50 60 70 
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Figure 5: Solution shape at t = 98 for the equation Uxxt + ^u^^Ux + UxxxU = 0. 




90 100 



Figure 6: Stable solution for the equation Uxxt+'iUxxUx+UxxxU — ^Uxxxx = 0, e = 0.01 
under the Gaussian initial condition. This figure is almost same as Figure ^ But, 
for large e they are quite different, and for negative e, the corresponding solution 
blows up. 
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Second, we give traveling wave solution for the 5th-order equation (|ip. Set 
u = ^ = X — ct {c is some constant speed to be determined), then after 

substituting this setting into equation (P we obtain 



12 

7 = y, c 



336 
625' 



So, the 5th-order equation (P has the following traveling wave solution 

, 336 ,12 



(69) 



(70) 



Although at each time the solution (|70D has singular point at x = —^t, this 5th- 



order PDE has the smooth and stable traveling wave solution under the Gaussian 
initial condition (see Figure p. 



equation u_t = -(u'^{-2/3})_{xxxxx} 




60 80 100 120 140 160 180 200 



Figure 7: This is the stable solution for the 5th-order equation Ut = d^u '^^^ under 
the Gaussian initial condition. 
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equation u_t = (u'^{-0.5})_{xxx} 




0.1 I 1 ' ' ' 1 1 1 1 1 1 

20 40 60 80 100 120 140 160 180 200 

X 



Figure 8: This is the stable solution for the HD equation Ut = d^u under the 
Gaussian initial condition. 



So, the figure ^ of the equation Ut = d^u has a slight difference from the 
figure H of the HD equation ut = d^u'^/"^. 

Third, we give traveling wave solution for the new integrable 7th-order equation 
@. Set u = ^ = X — ct {c is some constant speed to be determined), then we 
have 

_ 18 _ 31680 

^■y ""'Turn' ^ ^ 

So, the 7th-order equation has the following traveling wave solution 

31680 ,18 



Furthermore, we propose the following new equation: 

Ut = a^u""*/", / > 1, n ^ 0, m,neZ. (73) 
This equation has the following traveling wave solution 

u{x,t) = (x - ct)-"('-^)/('"+"\ (74) 

m T-r /mil 
c=-[[' 



k=l 



Apparently, if mn + > this solution has singularity at xq = cto (to is some 
time), and if mn + < this solution is a smooth traveling wave solution. 
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Remark 3 Here are the special cusp-like traveling wave solutions 

u{x,t) = {x--t)-^'^ (75) 
9 

and 

u{x,t) = {x+^^^tr^'^ (76) 
for the Harry-Dym equation Ut = d^{u~^^'^) and the 5th order equation Ut = d^iu""^^^). 



7 Conclusions 

In section 5, we obtain a parametric solution ( [55| ) of the 5th-order equation (|I]). This 
parametric solution can not include its traveling wave solution u = [x + |i|^)~~ 
because the constrained relation 

336 



does not hold, where 



{kqi{x,t),p3{x,t)) = + 



dlq. = -(- + i^^rA,: 



Traveling wave solutions u = {x + |||t) ^ for equation Ut = d^u and u = 
(x — n764°9 ^)~^ for equation Ut = d^ ^^ ^)xx-2{u 5)^ singular at some certain 
ponits X with the different time t. That is, this singularity travels with the time t. 
In the case of n{m + n) > 0, the traveling wave solution (|74D for general equation 
( [73D is also matching this case. 

A natural question arises here: is the equation ut = d\.u~^/^ integrable for all 
/ > 1, m, ri, G Z or for what kind of / > 1, m, G Z it is integrable? So far, for the 
cases: / = 3,m=l,n = 2 and / = 5, m = 2, n = 3, we know that it is integrable. 

The Harry-Dym equation has the cusp-like traveling wave solution u{x,t) = 
(x — |t)~^/^, but this is not cusp soliton which Wadati described this in Ref. | |22|| , 
because the traveling wave solution is singular, but the cusp is continuous. 

If we consider other constraints between the potential and the eigenfunctions, 
then we can still get the parametric solutions for other two equations (El) and IM). 
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